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Abstract: We combine Stein's method with Malhavin calculus in order to obtain explicit bounds in 
the multidimensional normal approximation (in the Wasserstein distance) of functionals of Gaussian 
fields. Our results generalize and refine the main findings by Peccati and Tudor (2005), Nualart 
and Ortiz-Latorre (2007), Peccati (2007) and Nourdin and Peccati (2007b, 2008); in particular, they 
apply to approximations by means of Gaussian vectors with an arbitrary, positive definite covariance 
matrix. Among several examples, we provide an application to a functional version of the Breuer- 
Major CLT for fields subordinated to a fractional Brownian motion. 
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Resume: Nous expliquons comment combiner la methode de Stein avec les outils du calcul de 
Malliavin pour majorer, de maniere explicite, la distance de Wasserstein entre une fonctionnellc d'un 
champs gaussien donnee et son approximation normale multidimensionnelle. Notre travail generalise 
et affine des resultats anterieurs prouves par Peccati et Tudor (2005), Nualart et Ortiz-Latorre 
(2007), Peccati (2007) et Nourdin et Peccati (2007b, 2008). Entre autres exemples, nous associons 
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du mouvement brownicn fractionnaire. 
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1 Introduction 

Let Z ~ o/K(0, 1) be a standard Gaussian random variable on some probability space 
(r2,J^, P), and let F be a real- valued functional of an infinite-dimensional Gaussian field. 
In the papers |12[ [T^ it is shown that one can combine Stein's method (see e.g. [3], [21] 
or [22]) with Malhavin calculus (see e.g. [H]), in order to deduce explicit (and, sometimes, 
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optimal) bounds for quantities of the type d(F,Z), where d stands for some distance be- 
tween the law of F and the law of Z (e.g., d can be the Kolmogorov or the Wasserstein 
distance). The aim of this paper is to extend the results of [12\ [TH] to the framework of 
the multidimensional Gaussian approximation in the Wasserstein distance. Once again, our 
techniques hinge upon the use of infinite-dimensional operators on Gaussian spaces (like 
the divergence operator or the Ornstein-Uhlenbeck generator) and upon an appropriate 
multidimensional version of Stein's method (in a form close to Chatterjee and Meckes |4j, 
but see also Reinert and Rollin (2^). As a result, we will obtain explicit bounds, both 
in terms of Malliavin derivatives and contraction operators, thus providing a substantial 
refinement of the main findings by Nualart and Ortiz-Latorre pLSj and Peccati and Tudor 
|18j . Note that an important part of our computations (see e.g. Lemma l3.7p are directly 
inspired by those contained in [I5J: we shall indeed stress that this last reference contains a 
fundamental methodological breakthrough, showing that one can deal with (possibly multi- 
dimensional) weak convergence on a Gaussian space, by means of Malliavin-type operators 
and "characterizing" differential equations. See [TT] for an application of these techniques to 
non-central limit theorems. Incidentally, observe that the paper [T7], which is mainly based 
on martingale-type techniques, also uses distances between probability measures (such as 
the Prokhorov distance) to deal with multidimensional Gaussian approximations on Wiener 
space, but without giving explicit bounds. 

The rationale behind Stein's method is better understood in dimension one. In this 
framework, the starting point is the following crucial result, proved e.g. in [21]. 

Lemma 1.1 (Stein's Lemma). A random variable Y is such that Y Z ~ ^(0, 1) if and 
only if, for every continuous and piecewise continuously differentiable function / : M ^ M 
such that E\f'{Z)\ < cxd, one has 

E[f'{Y) - YfiY)] = 0. (1.1) 

The fact that a random variable Y satisfying (jl.ip is necessarily Gaussian can be proved 
by several routes: for instance, by taking / to be a complex exponential, one can show that 
the characteristic function of y, say ip{t), is necessarily a solution to the differential equation 
V''(t) + tijjit) = 0, and therefore il}{t) = exp(— 1^/2); alternatively, one can set f{x) = x", 
n = 1,2, and observe that (|l.ip implies that, for every n, one must have E{Y^) = E{Z'^), 
where Z ~ o/K(0, 1) (note that the law of Z is determined by its moments). 

Heuristically, Lemma [l . 1 1 suggests that the distance d{Y, Z), between the law of a random 
variable Y and that of Z ~ ^(0,1), must be "small" whenever E[f'{Y) - Yf{Y)] ~ 0, 
for a sufficiently large class of functions /. In the seminal works [21^ I22j. Stein proved that 
this somewhat imprecise argument can be made rigorous by means of the use of differential 
equations. To see this, for a given function g : M ^ M, define the Stein equation associated 
with g as 

g{x) - E[g{Z)]=h'{x) - xh{x), Vx G R, (1.2) 

(we recall that Z ~ ^(0, 1)). A solution to ()1.2p is a function h which is Lebesgue-almost 
everywhere differentiable, and such that there exists a version of h' satisfying ()1.2p for every 
X E M. If one assumes that g G Lip(l) (that is, if Hfi'llLip < Ij where || • \\Lip stands for the 
usual Lipschitz seminorm), then a standard result (see e.g. [22]) yields that (II. 2p admits a 
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solution h such that ||/i'||oo < 1 and ||/i"||oo < 2. Now recah that the Wasserstein distance 
between the laws of two real-valued random variables Y and X is defined as 

d^{Y,X)= sup \E[g{Y)] - E[g{X)]\ , 

geLip(l) 

and introduce the notation = {/ '■ 1 1 /'I loo < 1 1 /"I loo < 2}. By taking expectations on 
the two sides of (|1.2p . one obtains finally that, for Z ~ ^(0, 1) and for a generic random 
variable Y, 

d^{Y,Z)< sup \E[f'{Y)-Yf{Y)]\, (1.3) 

thus giving a precise meaning to the heuristic argument sketched above (note that an anal- 
ogous conclusion can be obtained for other distances, such as the total variation distance 
or the Kolmogorov distance - see e.g. [5] for a discussion of this point). We stress that 
the topology induced by on probability measures on M, is stronger than the topology 
induced by weak convergence. 

The starting point of \\.2\ [T3] is that a relation such as (jl.3p can be very effectively 
combined with Malliavin calculus, whenever y is a centered regular functional of some 
infinite dimensional Gaussian field. To see this, denote by DY the Malliavin derivative of Y 
(observe that DY is a random element with values in some adequate Hilbert space i^), and 
write L to indicate the (infinite-dimensional) Ornstein-Uhlenbeck generator (see Section [2] 
below for precise definitions). One crucial relation proved in [12j, and then further exploited 
in [13], is the upper bound 

d^^{Y,Z) <E\l- {DY,-DL-^Y)^\. (1.4) 

As shown in [12], when specialized to the case of Y being equal to a multiple Wiener-Ito 
integral, relation (jl.4p yields bounds that are intimately related with the CLTs proved in 
|15j and [16) . See [13] for a characterization of the optimality of these bounds; see again [T2] 
for extensions to non-Gaussian approximations and for applications to the Breuer-Major 
CLT (stated and proved in [^) for functionals of a fractional Brownian motion. 

The principal contribution of the present paper (see e.g. the statement of Theorem 13.51 
below) consists in showing that a relation similar to (|1.4p continues to hold when Z is 
replaced by a d-dimensional {d > 2) Gaussian vector F = {Fi, ...,Fd) of smooth functionals 
of a Gaussian field, and is the Wasserstein distance between probability laws on M*^ 
(see Definition 13.11 below). Our results apply to Gaussian approximations by means of 
Gaussian vectors with arbitrary positive definite covariance matrices. The proofs rely on 
a multidimensional version of the Stein equation ()1.2p . that we combine with standard 
integration by parts formulae on an infinite-dimensional Gaussian space. Our approach 
bears some connections with the paper by Hsu [10], where the author proves an hybrid 
Stein/semimartingale characterization of Brownian motions on manifolds, via Malliavin- 
type operators. 

The paper is organized as follows. In Section [2] we provide some preliminaries on Malli- 
avin calculus. Section [3] contains our main results, concerning Gaussian approximations by 
means of vectors of Gaussian random variables with positive definite covariance matrices. 
Finally, Section [J] deals with two applications: (i) to a functional version of the Breuer-Major 
CLT (see |2]), and (ii) to Gaussian approximations of functionals of finite normal vectors, 
providing a generalization of a technical result proved by Chatterjee in [3]. 
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2 Preliminaries and notation 



In this section, we recall some basic elements of Malliavin calculus for Gaussian processes. 
The reader is referred to [13] for a complete discussion of this subject. Let X = {X{h), h G 
^} be an isonormal Gaussian process on a probability space (0,.^, P). This means that X 
is a centered Gaussian family indexed by the elements of an Hilbert space S), such that, for 
every pair h,g £ S) one has that E[X{h)X{g)] = {h,g)fj. 

We let LP'{X) be shorthand for the space L^(r2, cj(X), P). It is well known that every 
random variable F G L?'{X) admits the chaotic expansion F = E{F) + Y^^=i ^n{fn) where 
the deterministic kernels fn, n > 1, belong to i^®" and the convergence of the series holds 
in LP'{X). One sometimes uses the notation /o(/o) = ^W]- Iii the particular case where 

:= L'^{T, a, /i), with (T, A) a measurable space and fi is a u-finite measure without atoms, 
the random variable In{fn) coincides with the multiple Wiener-ltd integral (of order n) of 
/„ with respect to X (see [HI Section 1.1.2.]). 

Let / € ^®^, g G iD®^ and Q < r < p /\ q. We define the rth contraction f ®r 9 f and 
g as the element of ^®(P+'3-2r-) gjyen by 

oo 

f®rg--= ^ (/, (g) . . . (g) ej^)ijcsr (g) (5, Ci, (g) . . . (g) ej^)j5®r, 

where {cfc. A: > 1} is a complete orthonormal system in S). Note that f ®q g = f ® g; also, 
if p = (7, then f ®p g = (/, g) f^isp- Note that, in general, f ®r 9 '^^ not a symmetric element 
of ^'^(p+'J-^r). ^j-^g canonical symmetrization of / (g)r 5 is denoted by f®r9- We recall the 
product formula for multiple stochastic integrals: 

pAg / \ / \ _ 

WM = Y.%) Qlp+,-2r{f®rg). 

Now, let be the set of cylindrical functionals F of the form 

F = ^{X{hi),...,X{hn)), (2.1) 

where n > 1, hi G and the function ip G '^°°(M") is such that its partial derivatives have 
polynomial growth. The Malliavin derivative DF of a functional F of the form (j2.ip is the 
square integrable i^-valued random variable defined as 

n 

DF = ^dMX{hi), . . .,X{K))hi, (2.2) 

i=l 

where di^p denotes the ith. partial derivative of f. In particular, one has that DX{h) = h 
for every h in S^. By iteration, one can define the mth derivative D^F of F £ J/', which is 
an element of L^(i7; ij®™"), for m > 2. As usual B™'^ denotes the closure of with respect 
to the norm || • \\m,2 defined by the relation ||F||^^2 = ^i^^] + EI^i E[\\D'Ff^^,]. 

Note that every finite sum of Wiener-Ito integrals always belongs to D*"'^ (Vm > 1). The 
Malliavin derivative D satisfies the following chain rule formula: if (p : ^ M is in 
(defined as the set of continuously differentiable functions with bounded partial derivatives) 



4 



and if [Fi, . . . , Fn) is a random vector such that each component belongs to B^'^, then 
^{Fi, . . . , Fn) is itself an element of D^'^, and moreover 

n 

D^iFi,...,Fn) = Y,9MF'i,---,Fn)DFi. (2.3) 
1=1 

The divergence operator 6 is defined as the dual operator of D. Precisely, a random 
element u of (O; ij) belongs to the domain of 6 (denoted by DomJ) if there exists a constant 
Cu satisfying \E[{DF,u) < Cn||-F||L2(f2) for every F £ in this case, the divergence of n, 
written 8{u)^ is defined by the following duality property: 

E\Fb(u)\ = E\{DF, u)s^, VF G B^'^. (2.4) 

The crucial relation (j2.4|) is customarily called the (Malliavin) integration by parts formula. 

In what follows, we shall denote by T = {Tt : t > 0} the Ornstein-Uhlenheck semigroup. 
We recall that, for every t > and every F £ L'^{X), 

oo 

Tt{F) = Y,e-^'Jn{F), (2.5) 

n=0 

where, for every n > and for the rest of the paper, the symbol J„ denotes the projection op- 
erator onto the nth Wiener chaos, that is onto the closed linear subspace of L?'{X) generated 
by the random variables of the form Hn{X{h)) with h £ f) such that \\h\\fj = 1, and Hn the 
nth Hermite polynomial defined by ()4.ip . Note that T is indeed the semigroup associated 
with an infinite-dimensional stationary Gaussian process with values in R-^, having the law 
of X as an invariant distribution (see e.g. [14, Section 1.4] for a more detailed discussion of 
the Ornstein-Uhlenbeck semigroup in the context of Malliavin calculus; see Barbour [Ij for a 
version of Stein's method involving Ornstein-Uhlenbeck semigroups on infinite-dimensional 
spaces; see Gotze [9j for a version of Stein's method based on multi-dimensional Ornstein- 
Uhlenbeck semigroups). The inEnitesimal generator of the Ornstein-Uhlenbeck semigroup 
is noted L. A square integrable random variable F is in the domain of L (noted DomL) if 
F belongs to the domain of 5D (that is, if F is in D^'^ and DF £ Dom5) and, in this case, 
LF = —5DF. One can prove that LF is such that LF = — 'Yln^=Q^Jn{F). As an example, 
if F = Iq{fq), with fq £ ij®'?, then LF = -qF. Note that, for every F £ DomL, one 
has E{LF) = 0. The inverse of the operator L acts on zero-mean random variables 
F £ L'^{X) as L~^F = - lJn{F). In particular, for every q>l and every F = Iq{fq) 

with fq £ one has that L'^F = -\F. 

We conclude this section by recalling two important characterizations of the Ornstein- 
Uhlenbeck semigroup and its generator. 

i) Mehler's formula. Let F be an element of L?'{X), so that F can be represented as an 
application from into M. Then, an alternative representation (due to Mehler) of the 
action of the Ornstein-Uhlenbeck semigroup T (as defined in (|2.5|) ) on F, is the following: 

Tt{F) = E[F{e-'a + ^1 - e-^'X)] \a=x. t > 0, (2.6) 

where a designs a generic element of M^. See Nualart jTll Section 1.4.1] for more details on 
this and other characterizations of T. 
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ii) Differential characterization of L. Let F G L'^{X) have the form F = f{X{hi), ...,X{hd)), 
where / G '^^(R'^) has bounded first and second derivatives, and hi G $j, i = 1, ...,d. Then, 

i,j = l ' 1=1 * 

See Propositions 1.4.4 and 1.4.5 in [14J for a proof and some generahzations of (j2.7p . 

3 Stein's method and Gaussian vectors 

We start by giving a definition of the Wasserstein distance, as well as by introducing some 
useful norms over classes of real-valued matrices. 

Definition 3.1. (i) The W^asserstem distance between the laws of two M'^-valued random 
vectors X and Y, noted d-w{X,Y), is given by 

dw{X,Y):= sup \E[g{X)] - E[g{Y)]\, 

where indicates the class of Lipschitz functions, that is, the collection of all func- 
tions (7 : M'^ — > M such that ll^Hiip := sup — 1| 1| < oo (with || • \\-^d the usual 

Euclidian norm on 



(ii) The Hilbert- Schmidt inner product and the Hilhert- Schmidt norm on the class of d x d 
real matrices, denoted respectively by (•, ■)h.s. and || • \\h.s.i are defined as follows: for 
every pair of matrices A and {A,B)h.s. '■= Tr{AB^) and ||A||//.s'. := {A, A)h.s.- 

(iii) The operator norm of a dx d matrix A over M is given by ||j4||op := sup||^||^^=;^ || AxUjgd. 

Remark 3.2. 1. For every d > 1 the topology induced by dw, on the class of all probabil- 
ity measures on M'^, is strictly stronger than the topology induced by weak convergence 
(see e.g. Dudley ^ Chapter 11]). 

2. The reason why we focus on the Wasserstein distance is nested in the statement of the 
forthcoming Lemma 13.31 Indeed, according to relation (j3.4p . in order to control the 
second derivatives of the solution of the Stein equation (j3.3|) associated with g, one 
must use the fact that g is Lipschitz. 

3. According to the notation introduced in Definition 13. l( ii). relation ()2.7p can be rewrit- 
ten as 

LF = {C, Hessf{Z))H.s. - {Z, Vf(.Z))^,, (3.1) 

where Z = (X{hi), X{h(i)), and C = {C{i,j) : i,j = 1, d} is the d x d covariance 
matrix such that C{i,j) = E{X{hi)X{hj)) = {hi,hj)s^. 
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Given a d x d positive definite symmetric matrix C, we use the notation ,Ad{0,C) to 
indicate the law of a d-dimensional Gaussian vector with zero mean and covariance C. 
The foUowing result, which is basically known (see e.g. ^ or [20])) is the d-dimensional 
counterpart of Stein's Lemma ll.ll In what follows, we provide a new proof which is almost 
exclusively based on the use of Malliavin operators. 

Lemma 3.3. Fix an integer d > 2 and let C = {C{i,j) : i,j = 1, d} be a d x d positive 
definite symmetric real matrix. 

(i) Let Y be a random variable with values in W^. Then Y ~ .yV^{Q,C) if and only if, 
for every twice differentiable function / : M'^ — > M such that E\{C ^Yiess f {Y)) h^s] + 
E\{Y,V f{Y))^d\ < oo, it holds that 

E[{Y, V/(y))K. - {C, Ressf{Y))H.s.] = 0. (3.2) 

(ii) Consider a Gaussian random vector Z ~ o/^^(0, C). Let (7 : M'^ ^ M belong to '^^(R'^) 
with first and second bounded derivatives. Then, the function Uo{g) defined by 

Uog{x) := -E[g{^ftx + VT^tZ) - g{Z)]dt 

is a solution to the following differential equation (with unknown function f): 

g{x) - E[g{Z)] = {x, V/(x))r, - (C, Ressf{x))H.s., x G M"^. (3.3) 
Moreover, one has that 

sup ||Hessf7ofi'(x)||//.s. < ||C"^||op ||C||^/^ hWup- (3.4) 

Remark 3.4. 1. If C = fx^Id for some o" > (that is, if Z is composed of i.i.d. centered 
Gaussian random variables with common variance equal to o"^), then 

l|C~"^l|op \\C\\l^ = Wa^'^ldWop ll-T^Idllo^^ = fT~^ 

2. Unlike formulae (jl.ip and ()1.2p (associated with one-dimensional Gaussian approxima- 
tions) the relation p.2p and the Stein equation ()3.3p involve second-order differential 
operators. A discussion of this fact is detailed e.g. in [U Theorem 4]. 

Proof of Lemma \3.S\ We start by proving Point (ii). First observe that, without loss 
of generality, we can suppose that Z = (Zi,...,Z^) := {X{hi), ...X{hci)), where X is an 
isonormal Gaussian process over = W^, the kernels hi belong to ^ (z = l,...,(i), and 
{hi,hj)sj = E{X{hi)X{hj)) = E{ZiZj) = C{i,j). By using the change of variable 2u = 
— logi, one can rewrite Uog{x) as follows 

/•oo 

Uogix) = / {i?b(e-"x + ^1 - e-2-Z)] - E[g{Z)]}du. 
Jo 
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Now define g{Z) := g{Z) — E[g{Z)], and observe that g{Z) is by assumption a centered 
element of L^(X). For q>0, denote by Jq(g{Z)) the projection oig{Z) on tfie qth. Wiener 
chaos, so that jQ(g{Z)) = 0. According to Mehler's formula (j2.6|) . 



^[^(e-x + ^l-e-2"Z)]U=z - E[g{Z)] = S[5(e-"x + y^T^^Z)]\,=z = T^g{Z), 
where x denotes a generic element of W^. In view of (|2.5p . it follows that 

UogiZ) = / T^giZ)du = / V e-''^Jg{g{Z))du = V - J,(5(^)) = -L-'giZ). 

Since g belongs to ^^(M'^) with bounded first and second derivatives, it is easily seen that 
the same holds for Uog. By exploiting the differential representation (j3.ip . one deduces that 

{Z,VUog{Z))^, - {C,RessUog{Z))H.s. = -LUog{Z) = LL~^g{Z) = g{Z) - E[g{Z)]. 

Since the matrix C is positive definite, we infer that the support of the law of Z coincides 
with M"', and therefore (e.g. by a continuity argument) we obtain that 

{x,VUog{x))^d - {C,RessUog{x))H.s. = g{x) - E[g{Z)], 

for every x E W^. This yields that the function U^g solves the Stein's equation ()3.3p . 

To prove the estimate (|3.4p . we first recall that there exists a unique non-singular sym- 
metric matrix A such that = C, and that one has that A~^Z ~ cyfrf(0, I^)- Now write 
Uog{x) = h{A~^x), where 



h{x) 



/■I 1 

/ -E[gA{Vix + VT^A-'Z)-gA{A-'Z)]dt, 
Jo 



and gAix) = g{Ax). Note that, since A~^Z ~ ^/(^(O,!^), the function h solves the Stein's 
equation {x,Vh{x))^d — Ah{x) = gA{x) — E[gA{Y)], where Y ~ ^(0,1^). We can now use 
the same arguments as in the proof of Lemma 3 in [3] to deduce that 

sup ||HeSs/l(x)||//.S', < IbylllLip < ll^llopllffllLip- (3.5) 



On the other hand, by noting /ly^-i(x) = h{A~^x), one obtains by standard computations 
(recall that A is symmetric) that }iessUog{x) = Hess/i^-i(x) = A^^}iessh{A~^x)A~^, 
yielding 

sup ||Hess ?7o(7(x)||//,5. = sup \\A~^B.essh{A~^x)A~'^\\H.s. 

= sup \\A~^i{ess h{x)A~^\\H.s. 



< \\A ^ll^p sup ||Hess /i(x)||h.s'. (3-6) 

— 11^ Wop ll^llop ll^llLip (3-7) 

< IIC'^'^llop IIC'llop^ llsll-Lip- (3-8) 
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The chain of inequahties appearing in formulae (I3.6p - (l3.8p are mainly a consequence of the 
usual properties of the Hilbert-Schmidt and operator norms. Indeed, to prove inequality 
(|3.6|) we used the relations 

\\A~^Ress h{x)A~'^\\H.s. < W^'^Wop ||Hess/i(x)A"^||//.5'. 

< ||A"^||op ||Hess/i(x)||H.S'. ll^~"^l|op ; 

relation ()3.7p is a consequence of ()3.5p : finally, to show the inequality ()3.8p . one uses the 
fact that 

\\A-^\\op<^/lA^^^^^^7p = ^/^C^^P and \\A\\op < ^J\\AA\\op = ^J\\C\\op ■ 

We are now left with the proof of Point (i) in the statement. The fact that a vector 
Y ~ ^(0, C) necessarily verifies ()3.2p can be proved by standard integration by parts. 
On the other hand, suppose that Y verifies p.2p . Then, according to Point (ii), for every 
g G 'i^^(M'^) with bounded first and second derivatives, 

E{g{Y)) - E{g{Z)) = E{{Y,VUog{Y))^a - {C,RessUog{Y))H.s.) = 0, 

where Z ~ ^di^i C). Since the collection of all such functions g generates the Borel u-field 
on W^, this implies that Y Z, thus yielding the desired conclusion. □ 



The following statement is the main result of this paper. Its proof makes a crucial use 
of the integration by parts formula (|2.4p discussed in Section [2j 

Theorem 3.5. Fix d > 2 and let C = {C{i,j) : i,j = 1, ...,d} be a d x d positive definite 
matrix. Suppose that Z ~ ^d{0, C) and that F = {Fi, . . . , F^) is a -valued random vector 
such that E[Fi] = and Fi G B^'^ for every i = 1, . . . , d. Then, 



dy,{F,Z) < \\C-Xp\\C\\]l^^E\\C-^{DF)\\ls 



op 



\op 



\ 



^ E|(C(i,i) - {DF^-DL-^Fj}^)- 



(3.9) 
(3.10) 



where we write ^{DF) to indicate the matrix ^(DF) := {{DFi, —DL ^Fj)^ ■ ^ 1^ i, j 1^ d}. 

Proof. We start by proving that, for every g G '^^(M'^) with bounded first and second 
derivatives. 



\E[g{F)] - E[g{Z)]\ < WC^X^ \\C\\li^ \\g\\up^J E\\C - ^DF)rHs- 

To prove such a claim, observe that, according to Point (ii) in Lemma 13. 3| E[g{F)] 
E[g{Z)] = E[{F,VUog{F))^a - {C,}lessUog{F))H.s.]. Moreover, 

\E[{C,UessUog{F))H.s. - {F,VUog{F))^a]\ 

d d 
E C{i,j)df^Uog{F) - J2 mUog{F) 

i,j=l 1=1 
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E [Cii,j)df^UogiF)] -Y,E[{LL-'F,)dMF)] 

i,j=l i=l 
d d 

E [C{i,j)dlU^g{F)\ +YE[S{DL-'F,)diUog{F)] 

i,j=l i=l 
d d 

Y E [Cii,j)dfjUogiF)]-J2E[{DidiUo9{F)),-DL-'Fi)^ 

i,j=l i=l 
d d 

Y E [Cii,j)dfjUog{F)] -Y.^ [d]iUog{F){DF,,-DL-'F,)^ 

i,j=l i,j=l 



Y E [dlUog{F){C{i,j) - {DFi, -DL-'Fj)sj)] 

■id=i 

E{Ress Uog{F), C - ^{DF))h.s. \ 



since E{Fi) = 0) 
(since 5D = —L) 

(by (I23D) 
(by dZI 



< 
< 



E\\RessUQg{F)\\jjgJE\\C - ^{DF)\\\jg (by the Cauchy-Schwarz inequality) 



\C-Xv ||C||i/2 \\g\\L^p^E\\C-^{DF)\\l, (by m)- 



To prove the Wasserstein estimate ()3.9p . it is sufficient to observe that, for every globally 
Lipschitz function g such that H^Hlip ^ 1; there exists a family {g^ : e > 0} such that: 

(i) for each e > 0, the first and second derivatives of g^ are bounded; 

(ii) for each e > 0, one has that H^elliip < H^Hiip; 



(iii) as e ^ 0, \\ge - g\\oo i 0. 
For instance, we can choose ge{x) = E\^g{x + y/eN)\ with N ~ ^(0,1^) 



□ 



Observe that Theorem 13.51 generalizes relation (II. 4p (that was proved in [12^ Theorem 
3.1]). We now aim at applying Theorem 13.51 to vectors of multiple stochastic integrals. 

Corollary 3.6. Fix d > 2 and 1 < qi < ■ ■ ■ < qd- Consider a vector F := (Fi, . . . , F^) = 
(/g^(/i), . . . , Iq^ifd)) with fi e for any i = 1 . . . , d. Let Z ~ ./(^(O, C), with C positive 
definite. Then, 



dw{F,Z) < \\C-'\\op \\C\\ 



1/2 
op 



l<i,j<d 



C{i,j)--{DFi,DFj)s^ 



Proof. We have — L ^Fj = Fj so that the desired conclusion follows from ()3.10p . 



(3.11) 



□ 



When one applies Corollarv 13.61 in concrete situations (see e.g. Section [J] below), one can 
use the following result in order to evaluate the right-hand side of ()3.1ip . 
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Lemma 3.7. Let F = Ip{f) and G = Ig{g), with f € Sj^P and g € 53®'' (p,q>'^ )■ Let a he 
a real constant. If p = q, one has the estimate: 



E 



< ia-pl{f,g)sjs>p)" 



^a-^-{DF,DG)^J 

+ y - 1)'' I !) - 2r)! (11/ fWl^,^ + ||g g\\l^,^). 

r=l ^ 



On the other hand, if p < q, one has that 

1 V 

--{DF,DG)^j 



g-1 



Remark 3.8. 1. Recall that E{l,{f)I,{g)) = | ^^^^'^ otWise. 

2. In order to estimate the right-hand side of ()3.1ip . we see that it is sufficient to asses 
the quantity ®r fiW^^miqi-r) for any i G {1, . . . , d} and r G {1, . . . , — 1} on the 
one hand, and (/j, fj)f,<sii for any I < i, j < d such that qi = qj on the other hand. 

Proof of Lemma \3. ?| (see also [151 Lemma 2]). Without loss of generality, we can assume that 
= L'^{A,s^ , where {A,^/) is a measurable space, and ^ is a cr-finite and non-atomic 
measure. Thus, we can write 



{DF,DG)s, = pq{Ip-i{f),Ig^i{g))^=pq Ip-i{f{;t))l,^i{g{;t))f,{dt) 



pAq—1 



pq Yl 

-'^ r=0 

pAg— 1 

pq f^- 

r=0 



p-l 

r 

p-l\fq-l 
r I \ r 



1 



r 



Ip+g-2-2r{fi:,t)^rg{-,t))fl{dt) 



pq 



Ip+q-2~2r (f'^r+ig) 



It follows that 



E 



1 



a - -{DF,DG)i 



^'+p'E'r=i{r-iy-H^r-l)\r-yip+Q-'^^y-\\f^rgr 



(3.12) 



Sj0(p+q-2r) 



if p < g, 



. (a-p!(/,5V.)'+p'E?=i(^-l)!'rD'(2p-2r)!||/®,<7ll|«M P 
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If r < p < g then 

l|/'<5r-5ll^8)(p+9-2r) ^ 11/ "Xir 5'll^i8(p+9-2r) = (/ ®p-r /) ®g-r 
< \\f <^p-r f\\sj'S2r\\g (^q-r g\\sj^2r 

l{ r = p < q, then 

ll/^pS'll^^Cg-p) ^ 11/ ®p 5'll^®{<j-p) ^ ll/ll^®P ll^ 511^)82^. 

li r = p = q, then /®p5 = {f,g)sj<sp- By plugging these last expressions into ()3.12p . we 
deduce immediately the desired conclusion. □ 

Let us now recall the following result, which is a collection of some of the findings con- 
tained in the papers by Peccati and Tudor [TS] and Nualart and Ortiz-Latorre |15j . 



Theorem 3.9 (See [13 EH]). Fix d > 2 and let C = {C{i,j) : i,j = l,...,d} be a d x d 
positive definite matrix. Fix integers 1 < qi < ■ ■ ■ < qd- For any n > 1 and i = 1, . . . ,d, let 
/j^"^ belong to ij®"^*. Assume that 

= (Ff ), . . . , ) := (/,,(/(")), . . . , n > 1, 

is such that 

lim F[F/"^Ff ^] = C{i,j), l<i,j <d. (3.13) 
Then, as n —> oo, the following four assertions are equivalent: 

(i) For every I <i < d, f/"^ converges in distribution to a centered Gaussian random 
variable with variance C{i,i). 

(ii) For every 1 < i < d, E (f/"^ 3C{i,i)'^. 

(Hi) For every 1 < i < d and every I < r < qi — 1, H//"-* 0r fi^^\\^m(g—r) 0. 

(iv) The vector F^""^ converges in distribution to a d- dimensional Gaussian vector ^^{0, C) . 

Moreover, if C{i,j) = 5ij, where 6ij is the Kronecker symbol, then either one of conditions 
(i)~(iv) above is equivalent to the following: 

(v) For every 1 < i < d, ||Df/"^||^ q^. 

We conclude this section by pointing out the remarkable fact that, for vectors of mul- 
tiple Wiener-Ito integrals of arbitrary length, the Wasserstein distance metrizes the weak 
convergence towards a Gaussian vector with positive definite covariance. Note that the next 
statement also contains a generalization of the multidimensional results proved in to 
the case of an arbitrary covariance. 
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Proposition 3.10. Fix d > 2, let C be a positive definite d x d symmetric matrix, and let 
I < qi < ■ ■ ■ < Qd- Consider vectors 

F^"-^ := . . . ,Ff ) = . . . )), n > 1, 

with /f"^ G ^©'f^ for every i = 1 . . . ,d. Assume moreover that F^"^ satisfies condition iS.lS]) . 
Then, as n —> oo, the following three conditions are equivalent: 

(a) dwiF^'^lZ) 0. 

(b) For every I < i < d, q^^\\DFj^'^^\\j^ C{i,i) and, for every I < i j < d, 
{DFt\ -DL-^F^-^),, = qj\DFl-\DF\-^)^ ^ C{^,J). 

(c) -F^"-* converges in distribution to Z c/(^(0,C). 

Proof. Since convergence in the Wasserstein distance implies convergence in distribution, 
the implication (a) (c) is trivial. The implication (b) — > (a) is a consequence of relation 
(j3.1ip . Now assume that (c) is verified, that is, F^") converges in law to Z ~ ^(0, C) as n 
goes to infinity. By Theorem l3.9l we have that, for any i £ {1, . . . ,d} and r G {1, . . . , — 1}, 

By combining Corollarv 13.61 with Lemma 13.71 (see also Remark 13. 8r 2)). one therefore easily 
deduces that, since (I3.13P is in order, condition (b) must necessarily be satisfied. □ 

4 Applications 

4.1 Convergence of marginal distributions in the functional Breuer-Major 
CLT 

In this section, we use our main results in order to derive an explicit bound for the celebrated 
Breuer-Major CLT for fractional Brownian motion (fBm). We recall that a fBm B = {Bt : 
t > 0}, with Hurst index H £ (0, 1), is a centered Gaussian process, started from zero and 
with covariance function E[BsBt) = R{s,t), where 

R(s, t) = ^ {t^^ + -\t- s\^^) ; s, i > 0. 

For any choice of the Hurst parameter H £ (0, 1), the Gaussian space generated by B can 
be identified with an isonormal Gaussian process of the type X = {X{h) : h € S^}, where 
the real and separable Hilbert space Sj is defined as follows: (i) denote by £' the set of all 
M- valued step functions on [0,oo), (ii) define Sj as the Hilbert space obtained by closing 
with respect to the scalar product (ljo,t], l[o,s])^ = R{t,s). In particular, with such a 
notation, one has that Bt = X(1[q_^]). The reader is referred e.g. to [Tl] for more details on 
fBm, including crucial connections with fractional operators. We also define p{-) to be the 
covariance function associated with the stationary process x i— > B^+i — B^ (x S M), that is 

p(^) := + 1|2H ^ 1^ _ i|2// _ 2|^|2H) ^ ^|2^ _ ^1 |^|2//-2_ 

2 \x\^oo 
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Now fix an integer q > 2, assume that H < 1 — ^ and set 

[nt\-l 



Sn{t) — 1= Hq{Bk+l - Bk), t > 

awn ^-^ 

^ k=0 



where Hq is the gth Hermite polynomial defined as 



Hqix) 



(4.1) 



and where a = y^ql Ylrez {'"')■ According e.g. to the main results in [2] or [8], one has the 
following CLT: 



{Sn{t), t > 0} 



f.d.d. 



standard Brownian motion, 



where 'f.d.d.' indicates convergence in the sense of finite-dimensional distributions. To 
our knowledge, the following statement contains the first multidimensional bound for the 
Wasserstein distance ever proved for {Sn{t), t > 0}. 

Theorem 4.1. For any fixed d > 1 and = to < ti < . . . < td, there exists a constant c, 
(depending only on d, H and {tQ,ti, . . . ,td), and not on n) such that, for every n > 1: 



( -1 

n 2 



Sn{ti) Sn{ti—l) 
\/ti — ti—i 



;^d(0,ld) < cx { 



l<i<d 



n 



H-1 



2' 2q-2i 



Proof. Fix d > 1 and tQ = < ti < . . . < td- In the sequel, c will denote a constant 
independent of n, which can differ from one line to another. 
First, observe that 



MiLJ^t^ = i^^fM) with /f) 



1 



a^^Jti - ti-i 



lnti\-l 



In 



1}: 



fc=[nti_ij 

proof of Theorem 4.1, it is shown that, for any i G {1, . . . ,d} and r G {1, . . . ,qi 



Wfi"^ ®r /r^|L®2(,,-r) < C X < 



p{n) 



_ 1. 

n 2 



n 



H-l 



if ^ e ( 1 22^:31 



2' 2g-2J 



(4.2) 



Moreover, when 1 < i < j < d, we have: 

\,An) „(n)v I 
\{Ji ,Jj /fi®-?! 
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[ntij-l lntj\-l 
fc=[nii_ij Z=[ntj_iJ 



C 

n 

< c 



[nijj~[nti_ij-l 



^ [(K-J - 1 - r) A (Lnt,J - 1) - (K-_iJ - r) V (K_iJ)]p'?(r) 

r| = [ritj_ij-[ntij+l 

[ntjj - [ntj_ij - 1 



n 



^ 1^(^)1'/ = 0(^2gH-2,+i)^ asn^oo, (4.3) 

kl>L"*j-iJ-L"*d+i 
the last equality coming from 

^ |p(r)|^ = 0( E kl"""""^) = 0(iV"^^-'^+'), as iV ^ oo. 

|r|>Af \r\>N 

Finally, by combining (j4.2p . ()4.3p . Corollary 13.61 and Lemma |3.7^ we obtain the desired 
conclusion. □ 



4.2 Vector-valued functionals of finite Gaussian sequences 

Let Y = (Yi, ...,Yn) ~ o4^(0, I„), and let / : M*^ — > M be an absolutely continuous function 
such that / and its partial derivatives have subexponential growth at infinity. The following 
result has been proved by Chatterjee in [3], in the context of limit theorems for linear 
statistics of eigenvalues of random matrices. We use the notation d^y to indicate the total 
variation distance between laws of real valued random variables. 

Proposition 4.2 (Lemma 5.3 in [3]). Assume that the random variable W = f{Y) has zero 
mean and unit variance, and denote hy Z ^ ^(0, 1) a standard Gaussian random variable. 
Then d^iyiW, Z) < 2Var(T(y))^/^, where the function T(-) is defined as 

In what follows, we shall use Theorem 13.51 in order to deduce a multidimensional gener- 
alization of Proposition 14.21 (with the Wasserstein distance replacing total variation). 

Proposition 4.3. Let Y ~ ^n{0,K), where K = {K{i,l) : i,l = is a n x n 

positive definite matrix. Consider absolutely continuous functions fj : M" -^M, j = 1, ...,d. 
Assume that each random variable fj(Y) has zero mean, and also that each function fj 
and its partial derivatives have subexponential growth at infinity. Denote by Z ^ ^d{0, C) a 
Gaussian vector with values in Mf^ and with positive definite covariance matrix C = {C{a,b) : 
a, 6 = 1,..., 4. Finally, write W = {Wi, ...,Wd) = {fi(Y), fa{Y)). Then, 



dw{W,Z) < \\C-%p ||C||^/2 



Y,E[iC{a,b)-T,,iY)y] 
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where the functions Tab{-) are defined as 



Tab{y) = 




Proof. Without loss of generality, we can assume that Y = (Yi, 1^) = X{hn)), 
where X is an isonormal Gaussian process over some Hilbert space f), and {hi,hi)sj = 
E{X{hi)X{hi)) = K{i,l). According to Theorem 13.51 it is therefore sufficient to show that, 
for every a,b = \,...,d, TabiY) = {DWa, —DL~^Wb)f}. To prove this last claim, introduce 
the two ^-valued functions @a{y) and 0{,(y), defined for y S M"^ as follows: 

®-(y) = E§|(2^)^^ ®biy) = £ ^^p^[E[^{V'ty + Vi^tY)]h,'^dt. 

By using ()2.2p . it is easily seen that QaiY) = DWa- Moreover, by using e.g. formula 
(3.46) in [12], one deduces that eb{Y) = -DL-^Wb- Since Tab{Y) = (Ga(y), Gb(y))^, the 
conclusion is immediately obtained. □ 

By specializing the previous statement to the case n = d and fj{y) = yj, j = 1, ...,d, one 
obtains the following simple bound on the Wasserstein distance between Gaussian vectors 
of the same dimension (the proof is straightforward and omitted). 

Corollary 4.4. Let Y ~ o/f^(0,i^) and Z ./1^(0,C), where K and C are two positive 
definite covariance matrices. Then dYj{Y,Z) < Q{C,K) x ||C — ir||ff.5., where 

Q{C,K) :=min{||C-i||op \\C\\%\\\K'X,p \\K\\%^}. 
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